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1 INTRODUCTION 



Previous work 

iLvnden-Bell fc Boilvl 
effect of twist on a magnetic 
\ 19921 ) studied the case of a force-free magnetic dipole 
that is s heared at the stellar s urface , expands and opens 
partially. iLvnden- Bell fc Boilvl l| 19941 ) studied the case of 
a quadruple force-free magnetic field anchored on a heavy 
conductor that expands outwards because of differential ro- 
tation at the base and found that after finite twist the field 
opens tota lly and takes the topology of a split monopole. 
[Ahl (|1994| ) studied a 2-D model of arcade topology where 
the magnetic field expands because of the motion of the 
foot points at the surface where the magnetic field emerges. 

In this paper we study the problem of force free mag- 
netic arcades of azimuthal sym metry in three dimensions us - 
ing the methods developed by ILvnden-Bell fc Boilvl (|l994l ) 
and we seek for analytical solutions confined between two 
cones in order to form an arcade. It is not very likely that 
we can find solutions that will provide essential physical in- 
sight by just a numerical solution of the partial differential 
equation. It proves more interesting to design a problem that 
leads to analytical solutions. The method we are going to use 
can be described as a "backwards method", namely we are 
going to assume axial symmetry and self-similarity. Then 
we are going to evaluate how much differential rotation we 
must impose at the surface in order to achieve this structure 
and what is the external pressure environment that balances 
the magnetic pressure. This method is somewhat unnatural 
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ABSTRACT 

We study self-similar analytical solutions for force-free magnetic field in azimuthal 
symmetry and arcade topology. We assume the existence of a poloidal magnetic field, 
anchored on a heavy spherical conductor. The field is changed by shearing the foot 
points of the arcade due to differential rotation. This rotation gives rise to a toroidal 
component in the magnetic structure which reacts by expanding the poloidal flux 
outwards. This could be a slow process at the early stages, however it becomes very fast 
at the final stages when the poloidal flux expands to infinity. We address the question 
of the pressure environment confining the arcade, a pressure profile proportional to r~ 4 
is particularly interesting as it allows finite twist before the field expands to infinity. 
Finally, some time evolution estimates are made to demonstrate the limitations of this 
study. 

Key words: Sun: Magnetic Field; Stars: Magnetic Fields 



because it inverts the cause and the result, nevertheless it 
simplifies the problem and gives analytical solutions to a 
highly non-linear MHD problem. The alternative "forward 
method" of attacking the problem is more physical, it gives 
a wider class of solutions but the accuracy of its results is 
limited due to the forms of trial functions employed in the 
solution. In the "forward method" we take a given pressure 
environment, some flux emerging from the base and some 
differential rotation and we have to determine the shape of 
the magnetic cavity and the details of the magnetic field 
inside the cavity. 
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lAlvl l|l99 if ) has pointed out that there is a maximum 
energy that can be stored in a field whose foot points are 
fixed on a sphere and an upper bound is given by the field 
that goes straight outward to infinity. However in the case 
of a plane arcade in two dimensions no such bound exists 
and indeed such a two dimensional arcade can be sheared 
by any amount so that the field energy can be arbitrarily 
large. One of our aims is to demonstrate that arcades around 
the sun start as in the Aly 2-D case but once they have 
expanded enough to encounter the significantly larger areas 
available at great radii they rapidly increase in height and 
the flux expands to infinity after only a finite shear, unlike 
the 2-D case but in perfect conformity to Aly's bound. This 
upper energy bound has been proven for the case of a field 
occupying the space outside a sphere. Following Aly's proof 
step by step, it can be extended to the arcade case, where 
the field is confined by a sphere and two cones, provided that 
there is no normal component of the field into the bounding 
cones. Indeed the configuration contains a current sheet at 
the conical boundaries. The pressure is continuous, as the 
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The magnetic field is always divergence free, V ■ B = 
and we can express it in terms of the flux P(r,d), through 
the cup of a sphere which lies nearer to the axis than 9: 



B = - — VPxe^ + Vf (5) 

From iMestell (|l999f ) the toroidal component of the field 
can be expressed in the following way: 




Figure 1. The geometry of the problem. The field is confined in 
the space outside the sphere and between the two cones. 



external pressure is equal to the magnetic pressure at the 
boundary, for the arcade to be in equilibrium. 



2 THE PROBLEM 

In this section we state the physical properties of the configu- 
ration we are studying. A well defined framework is essential 
for analytical solutions. 



2.1 The Grad-Shafranov equation 

Consider a magnetic field of axial symmetry rising from and 
returning to the surface of a sphere but confined between 
two cones of semi-opening angles 9 1 and 9 2 (Fig. 1). The 
magnetic field is force-free and obeys the equations of ideal 
MHD. In what follows the symbols have their usual meaning. 
The force-free condition is 



Bj, = 



0(P) 
2nr sin 9 



(6) 



We substitute expressions (5) and (6) into equation (4) 
and rewrite the angular derivatives in terms of fi = cos 9. 
We write equation (4) in component form, therefore P(r,/j.) 
and /3(P) must obey the following relations 



, d 2 P , 



k = 



df3(P) 
dP ' 



(7) 



(8) 



The combination of (7) and (8) leads to the Grad- 
Shafranov equation in spherical coordinates under axial 



bnafranov equation in spnerical coordina tes ui 
svmmetrv l|Grad fc Rubinlll958l ; IShafranov|[l96r3 ) 



2 d 2 P^ n 2 ,d 2 P 



-r 2d -m P{ P). 



dP 



(9) 



Our objective is to solve equation (9) for P(r,jj,), for a 
certain form of the (3(P). 

2.2 Self-similar solutions 

The solution to the above partial differential equation de- 
pends on the choice of the form of /3(P). There is little hope 
to find an analytical solution to the above equation unless 
we determine j3{P). Nevertheless we can look for separable 
solutions which are self-similar in r. We express the flux P 
as a product function of a power law on r and a function 
determining the angular part f(fi). We normalise the maxi- 
mum of /(p) to unity using the total flux F; a is the radius 
of the sphere where the field emerges 



j x B = 0. 

The non-relativistic version of the induction equation, 
that is ignoring the displacement current, is 



_ _ 47T . 

V x B = — j. 

c J 

Therefore, substituting equation (2) to equation (1) 



(2) 



(V x B) x B = 0. 



(3) 



Equation (3) is equivalent to the condition that the curl 
of the magnetic field is parallel to the magnetic field, 



V x B = aB. 



(1) P(r,n) = F(-) f{fi). (10) 

The demand for self-similar solutions constrains to a 
power law dependence for the source function /3(P) on P 
jLvnden-Bell fc Boilvlll994h . We define f3(P) in the following 

way 

0(P) = Cl P\P\ 1/l . (11) 
Let Co = cii ?1/ ' i a; the angular part of P(r,fj,) satisfies 

1(1 + 1)/( M ) + (1 - M 2 )/"(M) = -cg(l + ])f(p)\f(p)\ 2/l - (12) 

Therefore the non-linear partial differential equation (9) 
(4) reduces to an ordinary differential equation, which is still 
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non-linear, subject to boundary conditions and the con- 
straint f(fi) 1. In the next sections we are going to in- 
vestigate solutions with arcade topology, but before that we 
will discuss the form of the fields that appear in these con- 
figurations. 

2.3 The fields 

According to the relations introduced in sections 2.1 and 
2.2, it is sufficient to determine the form of f(fi), by solv- 
ing equation (12). We then obtain detailed expressions for 
the components of the magnetic fields and consequently the 
structure of the field lines. The poloidal component of the 
field is the first term of the right hand side of expression (5), 
the toroidal is given by expression (6), and by virtue of (10) 
and (11) the fields become: 



13 r 



Bo = 



dm 



dfx 



2tt 

Fa l lr- 1 - 2 
2tt(1 - m 2 ) 1/2 

c Fa l r-'- 2 



/(/*), 



dr 
13, 



2tt(1- M 2 ) 1 /2- 
The field lines are determined by the relation: 

rd8 r : 



(13) 
(14) 
(15) 

(16) 



In order to evaluate the twist of a field line, we are 
using the equation involving the 6 and <f> components of the 
magnetic field 



(17) 



rdO r sm f 

W = ^f\f\ 1/r 

Thus we can integrate and evaluate the total twist on 
a field line that emerges from a point at r — a, n — /ii and 
ends at r — a and (j, = /X2. 



<f> 



c o I /I 



l/i 



I 



/' 



(18) 



Thus, if we know the form of f(fi) we have all the in- 
formation about the field lines. 



3 ARCADE SOLUTIONS 

In this section we study the magnetic arcades in two distinct 
cases. The first one is a single arcade bounded between two 
cones. We increase the twist of the field lines at the base 
and consequently the arcade expands radially. We evaluate 
the twist of the field lines and the energy that is stored in 
this configuration. The second one is a triple arcade. The 
boundary condition confines the edges of the external struc- 
ture, while the central arcade, which is initially identical to 
the single arcade we study in the first paradigm, has limited 
mobility. We are interested in this configuration in order to 
mitigate the effect of the artificial conical boundaries and 
have a comparison with the single arcade. 




Figure 2. The angular function f(fi) of the single arcade for 
I 9, 1, 0.2, 0.005, dashed-dotted, dotted, dashed and solid lines 
respectively. 



3.1 The single arcade 

We solve equation (12) to construct a model for a single ar- 
cade. The field emerges from a spherical surface of radius a 
and is constrained between two co-axial cones, whose ver- 
tices coincide with the centre of the sphere. Thus the bound- 
ary conditions of this configuration, applicable to equation 
(12) are f(fJ.i) = /(/te) = 0, where pi and \i2 are the cosines 
of the semi-opening angles of the cones. The angular func- 
tion f(j-i) is normalised to unity, thus we constrain the value 
of Co for given I. We start from I = 9 and Co = 0, then we 
solve the differential equation for I decreasing, so that the 
structure expands (Fig. 2). We evaluate the total twist of 
the structure, and we can estimate the differential rotation 
we need to impose at each field line in order to achieve this 
configuration (Fig. 3). 

For Co = the field is purely poloidal and the solution 
takes the following analytical form 



h = 



2.188 



(19) 



Where -Pg(/i) is the derivative with respect to /1 of the 
ninth Legendre polynomial. 

The choice of I — 9 in the initial case yields that the 
boundaries of the arcade are Hi = 0.478 and ^2 = 0.738. The 
other cases are solved numerically and the results appear at 
table 1. Whe n I — > we can use the a pproximate method 
developed bv lLvnden- Bell fc Boilvl (|1994| ) in order to obtain 
an analytical expression for f(fi) 



m = 1 - 



/ 



1+/ 



■ In cosh 



(j + l)cp 
,2 ^1/2 



(1 ^maxo) 1 



(fl — flmaxo) I • (20) 



The quantities appearing equation (20) are determined 
by the asymptotic behaviour of the differential equation 
(12). As cosh is an even parity function, in order to sat- 
isfy the boundary condition at fii and Umax obeys 
Mi - /immo = — (M2 — fJ-max) therefore ^i ma x0 = {fii + M2)/2. 
Then it is possible to find Co using again the boundary condi- 

^/ ; >), in 



tion f(ui) = 0, thus c = 2 "'"^" -cosfe- 1 (e (1 
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Horizontal Shear 
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Figure 3. The dependence of the horizontal shear on the power 
law index I. Only finite differential rotation is needed for the ar- 
cade to reach infinity with / = 0. 




Figure 4. The angular function f(fJ.) of the triple arcade with l- 
9, 1, 0.2, 0.1, solid, dashed, dotted, dashcd-dottcd respectively.. 
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/'(Ml) 


$ 


Energy 


9 


0.00 


0.613 


11.35 


0.00 


1.00 


4 


7.87 


0.613 


10.87 


-0.71 


1.97 


2 


8.22 


0.613 


10.35 


-1.30 


3.22 


1 


7.87 


0.612 


9.72 


-1.99 


4.68 


0.5 


7.35 


0.611 


9.08 


-2.79 


6.17 


0.2 


6.76 


0.610 


8.40 


-3.31 


7.18 


0.1 


6.47 


0.609 


8.08 


-3.61 


7.65 


0.0 


6.11 


0.608 


7.69 


-3.96 


8.12 



Table 1. The results of the numerical solution for equation (12). 
<E> is total twist about the axis of the star of the field lines emerging 
and ending at the boundaries of the arcade corresponding to P = 
and the energy is scaled to the energy of the purely poloidal 
arcade. Expansion to infinity occurs when the twist is somewhat 
over half a turn. 



the limit of Z — + it becomes Co = 2(1— ft^o) 1 ' 2 / (/X2 — /w). 
The boundary condition /(^a) = is satisfied because of 
the symmetry of the function / with respect to /J, m ax0- The 
numerical values for Co and fi m ax are on tables 1 and 2. 
A higher order approximation can be d one by applying the 
more sophisticated method developed bv lLvnden-Belll (l2006t ) 
(Appendix A), that computes the variation of Umax as it 
reaches its asymptotic value. The twist for the limiting case 
I — > can be evaluated analytically as described in the ap- 
pendix. 

The numerical solution of equation (12) reveals the 
properties of the magnetic configuration. The increase of the 
twist amplifies the energy stored in the arcade by two means, 
the ij> component of the magnetic field becomes stronger and 
the arcaded expands. As opposed to the plane-parallel model 
developed by Aly, finite twist leads to disconnection of the 
arcade into an outgoing and an incoming part. The B<f, com- 
ponent concentrates at the cone separating the incoming and 
the outgoing part of the arcade and a sheet current forms 
in the e$ direction. The energy stored in the configuration 
increases, but is always finite. The foot points of the field 
lines move weakly in latitude, but somewhat more signifi- 
cantly close to the centre of the arcade. 
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/'(w) 




Energy 
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0.00 


0.478 


0.738 


0.613 


11.35 


0.00 


1.00 


4 


7.88 


0.472 


0.733 


0.607 


10.89 


-0.71 


1.98 


2 


8.29 


0.465 


0.726 


0.600 


10.33 


-1.28 


3.20 


1 


8.03 


0.456 


0.716 


0.589 


9.74 


-1.94 


4.66 


0.5 


7.81 


0.445 


0.704 


0.577 


9.15 


-2.64 


5.88 


0.2 


7.16 


0.432 


0.688 


0.562 


8.55 


-3.17 


7.26 


0.1 


6.93 


0.426 


0.680 


0.554 


8.28 


-3.43 


7.79 


0.0 


6.69 


0.417 


0.668 


0.543 


7.94 


-3.74 


8.36 



Table 2. The numerical solution of the triple arcade. The energy 
refers to the energy of the central arcade and $ to the twist of the 
field line that lies in the central arcade and has P = 0. There is 
minor deviation from the results for the single arcade appearing 
at table 1 because the foot points now move a little in lattitudc. 

3.2 The triple arcade 

In order to mitigate the effect of the external boundaries on 
the arcade we consider three neighbouring arcades, where 
the central one has initially exactly the same structure as 
the arcade of the previous section. At the untwisted stage, 
that is I — 9 and Co = 0; the other two are determined 
by the Legendre Polynomial we chose. Then, following the 
same process the foot points are differentially rotated in the 
azimuthal direction leading to the expansion of the field. The 
advantage of this method is that the artificial boundaries 
have been removed from the arcade we are interested in and 
we have allowed limited mobility. The whole structure is 
confined again between two cones of cosines of semi-opening 
angles /13 and /m that determine the boundary condition 
/(0.165) = /(0.920) = 0, these correspond to the first and 
the fourth root of the expression of (1 — pi 2 )Pg(/i); we reserve 
Hi and /12 for the boundaries of the central arcade (table 
1). We have chosen again I — 9 in the initial case so that 
the solutions have three humps, and in order to be able to 
compare directly the triple to the single arcade (Fig. 4). 

Following the same process we gradually decrease I 
down to zero and we solve numerically the problem under 
the normalisation f m ux = 1. The results appear in table 2. 

The overall image is similar to the previous case, the 
arcade expands outwards as we increase the twist of the foot 



Self-Similar Magnetic Arcades 5 



+ Energy of the Poloidal Field 
+ Energy of the Toroidal Field 
Total Energy 



7- 
6- 

CD 
C 

3 

2- 

V- 




Figure 5. The energy of the the arcade. Initially, the energy is 
only due to the poloidal component. As the shear increases there 
is energy in the toroidal component as well, however, after some 
point the toroidal component moves to greater distances from the 
base and its contribution to the total energy decreases. In the final 
state of the split monopole the energy is only due to the poloidal 
component, as expected. 



Figure 7. The magnetic field lines of the triple arcade as seen 
from the side for I = 2, the twist has caused a significant expan- 
sion. 





Figure 6. The magnetic field lines of the triple arcade as seen 
from the side in the untwisted stage for 1 = 9. 



points and energy is stored in the magnetic configuration 
(Fig. 5). The angular size of the central arcade varies very 
weakly with I, although both foot points move towards the 
equator and the twist is concentrated in the centre of the 
arcade (Figures. 6-10). 

The solutions of the single and the triple arcade give 
the essential understanding to this problem. We see little 
deviation between the two models, therefore mobility in the 
boundaries does not alter significantly the general proper- 



Figure 8. The magnetic field lines of the triple arcade as seen by 
the side for I = 0.4, a little before the field lines open. 



ties, apart from some minor changes in the total energy and 
the total twist of the configuration. Therefore in future and 
more detailed studies we can neglect minor movements in 
latitude at the base and keep the boundaries fixed. 
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Horizontal Shear, R = 



Figure 9. The twist of the magnetic field lines of a surface of 
constant flux, for 1 = 2. 




Figure 11. The size of the arcade as a function of the horizontal 
shear in the 2 D and 3 D models. We start with two arcades of 
the same size, when the 3 D arcade becomes comparable to the 
size of the star their evolution changes dramatically. 



the same size in 2-D and we evaluate the expansion as a func- 
tion of the horizontal shearing (Figure 11). Initially, both the 
2-D and 3-D arcades, they coincide, however when the size 
of the 3-D arcade becomes comparable to that of the star 
their evolution changes drastically, the 2-D model continues 
on a linear increase of the size of the arcade with the hori- 
zontal shearing, whereas the expansion of the 3-D is much 
more rapid. Besides that, there is need for infinite amount 
of energy in order to force the plane parallel arcade to open 
and reach infinity as opposed to the finite energy demands 
of the 3-D arcade. 




Figure 10. The twist of the magnetic field lines for / = 0.4, the 
expansion of the configuration and the increase in the twist are 
obvious, compared to figure 9. 



3.3 Spherical and plane parallel models 

An important conclusion of the previous section is that after 
finite shear the arcade expands and releases its energy. It is 
essential to have a seed poloidal arcade, where we can am- 
plify its energy by almost a factor of te n by t wistin g before 
it erupts. This is in agreement to what lAlvl (1991) proved, 
that only finite energy can be stored in a 3-D force-free mag- 
netic structure. The question of 2-D plane parallel structures 
and 3-D axially s ymmetric arcades is discussed in general by 
lUzdenskvl (|2002h where the physical deviations of the two 
configurations are stressed. 

A similar version of th i s prob lem has been solved by 
iLow fc Wolfsonl ljl988h : [/Uvl ij 19941 ) where solar flares were 
described through a 2-D plane parallel model. In order to 
have a direct comparison we construct a poloidal arcade of 



4 THE PRESSURE 

In this section we address the question of the pressure envi- 
ronment of the arcade. We estimate the dependence of the 
magnetic pressure on r and how it varies with I. Afterwards 
we comment on the evolution in arbitrary environments, we 
do not deduce detailed solutions for the structure of the 
field inside the magnetic cavity, because an arbitrary envi- 
ronment does not allow analytical self-similar solutions and 
a detailed numerical model is beyond the scope of this paper; 
however based on physical arguments about pressure equi- 
librium at the edge of the magnetic cavity, we gain some 
insight about its shape. 

There is a physical reason that prevents the existence of 
self-similar solutions in the case of an arbitrary environment. 
If we include an arbitrary pressure we increase the external 
parameters and we are able to construct by the means of 
dimensional analysis a typical length scale; systems having 
a typical length scale do not allow self-similar solutions. 



4.1 Equilibrium with an external environment 

In order to confine this structure, it has to balance with an 
external pressure at its boundaries. The magnetic pressure 
Pm at /ii and [12 is only due to the B r component of the 
magnetic field as the other two components vanish. 
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and by setting p m ,o = ^ (l£ir f' i the pressure is 

Pm(r) = p m ,or~ 2l ~ 4 . (22) 

If the pressure of the gas surrounding the structure is 
obeying a power law p g oc , then in order to have equi- 
librium the indexes must satisfy 

C = 21 + 4 (23) 

It is rather unphysical to expect a natural environment 
to obey this varying pressure law. 

It makes better sense to study this structure subject to 
a given pressure. But we have to sacrifice the advantage of 
self-similar analytical solutions and seek a numerical solu- 
tion. It could be studied using m ethods similar to the ones 
developed by iLynden-Bell l|2006l) . 

4.2 Evolution under given pressure environment 

In section 4.1 we found that the magnetic pressure is pro- 
portional to r~ 2l ~ 4 , and in our arcade models we set I to de- 
crease from 9 down to 0. Consider an environment where the 
pressure is proportional to and it is in equilibrium with 
the arcade at the surface of the star, thusp s (a) = p m (a). For 
any distance n that is n > a the tendency of the arcade 
to expand or to contract in the angular sense depends on 
the sign of Q — 21 + 4 — £. If Q is positive the gas pressure 
will be greater than the magnetic pressure, therefore, the 
arcade is more compressed and may indeed be pushed down 
and smothered. If Q is negative the arcade expands in the 
angular direction, whereas for Q = 0, the opening angle of 
the arcade will remain constant and the self-similar solution 
will satisfy the boundary conditions. 

In the cases studied, we assumed that I varies from 
to 9, therefore it is interesting to describe the behaviour of 
an arcade that evolves in an environment of certain £ par- 
ticularly when C > 4 and £ = 4, the case of £ < 4 is not 
discussed here because the system does not pass through a 
stage of self-similar solutions. The choice of I = 9 initially 
constrains £ < 13, otherwise it would be impossible to con- 
fine the arcade at any stage. 

Let £ > 4, initially I = 9 and Q > therefore the gas 
pressure collimates the arcade. However, as I decreases, there 
is some value of I for which Q = 0, at this point there is no 
deviation from the self-similar solutions. Just after that Q 
becomes negative, the opening angle of the arcade increases 
and the arcade opens to infinity. Therefore the system re- 
leases its energy faster than it does in the fully self-similar 
model. 

For £ = 4, Q is always positive, unless I = 4, therefore 
the arcade is always confined by its boundaries and its open- 
ing angle decreases as we move away from the surface of the 
star, the final part of its evolution is the split monopole, at 
this stage the arcade obeys the solution for I — and it 
consists of two streams of field lines, one outgoing and the 
other incoming that connect at infinity. 



The process can be summarised in the following stages. 
We start with a magnetic structure where the gas pressure 
is stronger than the corresponding pressure of a self-similar 
arcade. We inject energy to the system by shearing the foot 
points, and the arcade expands so that it reaches a stage 
where self-similar solutions are permitted, at this stage, the 
arcade is confined between two cones and is described by 
the solutions of the previous section. If we shear a bit more 
the foot points the arcade deviates from the self-similar so- 
lution and in addition to the radial expansion it opens in 
the angular direction. The amount of shearing and energy 
needed to achieve that, is defined by the dependence of the 
pressure to the distance and can be roughly approximated 
by the energy of the self-similar arcades appearing on table 
1. 



5 TIME EVOLUTION 

The discussion of the problem did not include any time de- 
pendence. We are going to study this aspect of the problem 
in more detail. Using the equilibrium solutions we find the 
dependence of the differential rotation at the surface $ on I. 
We are going to assume that $ increases linearly with time, 
and then we are going to see how fast the arcades expand 
due to the displacement at the surface. At this point we are 
going to estimate the relation between twist and expansion 
for a given field line that corresponds to P = 1 /2P max . Let 
the angular displacement of the foot points of a field line 
determined by P be $p, and let the displacement imposed 
at the base increase linearly with time for this field line 

$ P = Q P t (24) 

A motion at the origin 5^ = asin8S(j> corresponds to 
an expansion Sr. The expansion velocity cannot exceed the 
Alfven velocity in the medium hosting the arcade. This 
is because the perturbation responsible for the expansion 
cannot travel faster than a magnetic wave. In order to 
make an evaluation of the velocity we assume that the field 
emerges from a spherical star of the size of the sun, its 
strength is 10 2 Gauss and the particle density of the plasma 
inside the arcade is 10 9 c m" 3 , corresponding to an Alfven 
velocity of 6 x 10 8 cm/ sec l|Priest fc Forbes!l2002l ). that is 2% 
of the speed of light. The angular velocity of the differential 
rotation is 10~ 7 rad/ sec. Then we compute the expansion 
of the arcade corresponding to this shearing and we expect 
that both take place in the same time scale. Therefore 
we set a limit for the validity of the model of quasi static 
solutions. In table 3 we present the corresponding velocities 
and the relation between the shear and the expansion. We 
conclude that under the above assumption the magnetic 
configuration cannot be studied with this method when 
the power law index becomes less than 0.1. At the final 
stages the expansion velocity would be extremely fast and 
the field would not have time to come to an equilibrium. 
The details of this critical point depend on the choice of 
the parameters of the problem. In any case no matter how 
slow the shearing is, the expansion at some point becomes 
infinitely fast, a choice of a very small £1 may only postpone 
this for later stages. Therefore the assumption of quasi 
static solutions is limited. This problem demonstrates that 
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t(xlO b sec) 


$ 


I 




v exp (cm/sec) 








9 


1.08 


0.0 


1.7 


0.52 


4 


1.19 


4.4 x 10 4 


3.3 


1.00 


2 


1.41 


9.7 x 10 4 


5.5 


1.66 


1 


2.00 


1.8 x 10 5 


8.1 


2.42 


0.5 


4.00 


5.4 x 10 5 


10.9 


3.28 


0.2 


32 


6.8 x 10 6 


12.0 


3.60 


0.1 


10 3 


6.4 x 10 s 


12.6 


3.77 


0.05 


10 6 





Table 3. The time evolution of the field line corresponding to 
P = Pl/2- 



expansion should take place in the same time scale. This 
time scale has to be shorter than the Alfven crossing-time, 
otherwise the message of the change of the shear at the base 
will not be transmitted to the top of the arcade. Therefore 
making a rough estimate we can set an upper limit to the va- 
lidity of this method, relating the shear velocity at the origin 
and other physical quantities of the magnetic structure. We 
concluded that there is a limit related to the physical prop- 
erties of the system, the static solutions cannot be valid for 
ever, therefore a dynamical study of the problem is essential. 



it is essential to seek for other solutions that include time 
evolution explicitly and create a fully dynamical model. 
This task seems far from the capacity of the current ana- 
lytical model, however there are 3-D dynamical simulati ons 
liMurrav. Hood.Moreno-Insertis. Galsgaard fe Archontisl 
I2006J ) giving some insight. 



6 CONCLUSIONS 

In this study we looked in detail the possibility of con- 
structing magnetic a r cades. The Grad -Shafranov equation 
l|Grad &; Rubinlll958l ; IShafranovl Il966:) accepts self-similar 
solutions confined between two cones and of arcade topol- 
ogy. In order to avoid numerical computations we make 
certain assumptions about the distribution of the toroidal 
component of the magnetic field inside the structure. The 
starting point is a poloidal current-free magnetic field 
whereas the endpoint is the split monopole. At the split 
monopole stage a sheet current forms at the boundary be- 
tween incoming and outgoing magnetic field, reconnection 
effects must be very important at t his stage of evolution 
l|Amari. Luciani. Aly fc Taggerl Il996h . A similar case stud- 
ied is that of the triple arcade, in this scenario we study a 
triple arcade focusing on the central arcade. The boundaries 
of the central arcade are allowed to move in latitude on the 
surface of the star. We conclude that there is no great differ- 
ence between the two models, thus in future work, focusing 
on the single arcade is sufficient. 

The comparison of the 2-D and 3-D models is a very 
interesting topic. Two initially similar arcades, one in 2-D 
plane parallel model and the second in 3-D azimuthal sym- 
metry, diverge as their size increases. The critical point for 
the divergence occurs when the 3-D arcade becomes compa- 
rable with the size of the star. 

We also address the question of the pressure. The an- 
alytical idealised solutions we obtain may only appear in 
certain pressure environments. If we look for a more real- 
istic pressure environment we have to allow the boundaries 
of the arcade to vary with altitude and not to be confined 
between the two cones, however we have to sacrifice self- 
similarity and make a fully numerical model that is beyond 
the scope of this paper. However the use of our solutions 
provides essential insight about the evolution of such struc- 
tures. 

Finally we studied the time evolution of this structure. 
We considered a field line corresponding to P\/2 of the flux 
and we studied its evolution. Assuming that we make this 
shear under constant velocity at the base we expect that the 
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APPENDIX: THE TOTAL TWIST 

The twist of a field line emerging at and sinking to lat- 
itudes corresponding to fii and \xi respectively is given by 
the integral of the equation (18). When I — > the angu- 
lar function /(/x) can be approximated by the relation (20), 
section 3.1. Let 



w = In cosh 



( (I + l)co , A 

"l 1 H-maxl)) ' 



f =0---nr[ w ) ~exp(-w), 

therefore the integrand of equation (18) becomes 
cof 1/l 



CO 



(1) 

(2) 
(3) 



Kl ^ Kl - M 2 ) cosh ( K J£l2) l/2 ifl ~ Mma:c0) ) 

The variation of the integrand of equation (18) due to 
f 1 (fi) is much stronger than the variation due to 1 — /i 2 , 
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since well away from fj, m axO f is zero. We put then 1 — fi 2 = 
1 — cLiO an d we integrate. 



This integration can be done analytically and gives 
2 

$ = =t- — — {arctanexp[ci(//2 - ^maio)] 

yji. — fJ-maxO(l + l) 

- arctanexp[ci(/ii - Umaxo)]}, (5) 

where ci = (I + l)c /(Vl = V-maxo)- When / — » 0, ci 
becomes pretty big, therefore the first term in the bracket is 
7r/2, whereas the second term is 0. Thus, the total twist is: 




$ = _ 



(6) 



'maxO 
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